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It is well known that pulse-like solutions of the cubic complex Ginzburg-Landau equation are un- 
stable but can be stabilised by the addition of quintic terms. In this paper we explore an alternative 
mechanism where the role of the stabilising agent is played by the parametric driver. Our analysis is 
based on the numerical continuation of solutions in one of the parameters of the Ginzburg-Landau 
equation (the diffusion coefficient c), starting from the nonlinear Schrodinger limit (for which c = 0). 
The continuation generates, recursively, a sequence of coexisting stable solutions with increasing 
number of humps. The sequence "converges" to a long pulse which can be interpreted as a bound 
state of two fronts with opposite polarities. 



I. INTRODUCTION 



A variety of noncquilibrium phenomena such as open flows in hydrodynamics, thermal convection in pure fluids 
and binary mixtures, processes in lasers and nucleation during the first-order phase transitions, can be modelled by 
the complex Ginzburg-Landau equations. (For review and reference, see e.g. 0, y, Li 0, Q-) Of primary importance 
are pulse-like solitary wave solutions, which represent localised structures widely observed in nonequilibrium systems. 
It is well known that in the cubic Ginzburg-Landau equation, 

ivb t + i-ii> + (1 - ic)^ xx + (2 - ig)\ip\ 2 ip = 0, (1) 

solitary waves are unstable for all positive c, 7 and real g. In the g > case, however, they can be stabilised by the 
addition of quintic terms: 

iipt + ijijj + (1 - ic)ip xx + (2 - ig)\ijj\ 2 i) = -{q r + iqi)\ip\ 4 ip, (2) 

with positive qi • This example suggests that solitary pulses can be stable in a more general class of Ginzburg- 
Landau equations where the zero solution undergoes a subcritical (rather than supercritical) bifurcation to a flat 
nonzero solution 0. In Eq.®, the terms with 7 and c account for linear homogeneous and nonhomogeneous losses, 
respectively, while the terms with g and qi describe the cubic gain and quintic dissipation. 

The present work deals with another equation of the Ginzburg-Landau type exhibiting the subcritical bifurcation, 
viz. the parametrically driven Ginzburg-Landau: 

iik + Hi> + (1 - *c)Vw + (2 - ig) IV-I 2 rl> = H*e 2iut . (3) 

Here, as in the positive c and 7 are the homogeneous loss and diffusion coefficients, respectively. The term 

with an asterisk (indicating the complex conjugation) represents the parametric driver. (The driver's amplitude h 
can be chosen positive without loss of generality.) When c = g — 0, Eq.© gives the parametrically driven damped 
nonlinear Schrodinger equation. This special case has been studied extensively; in particular, stable solitary waves 
Q and their stable complexes were found, and their bifurcations and supercritical dynamics analysed ^(|. The 
objective of the present work is to advance beyond the nonlinear Schrodinger limit. We will still keep 5 = but 
allow for a nonzero diffusion coefficient, c. As we will see, even such a minimal generalisation gives rise to a new 
phenomenology of localised solutions which includes the multistability of pulses and pulse- front transitions. 

Like the Ginzburg-Landau equations with intrinsic gain, the parametrically driven equation © arises in a wide 
ran ge of p hysical applications. These include nonlinear Faraday resonance in vertically vibrated layers of water 
|lllll2llT3| and nonlinear lattices [l4j |: commensurate- incommensurate transitions in convective systems 15]; waves in 
nematic and cholesteric liquid crystals in rotating magnetic fields |16| ; magnetisation waves in easy-plane ferromagnet 
exposed to microwave fields || ; domain walls in the easy-axis ferromagnet near the Curie temperature |17| and in the 
easy-plane magnet in the stationary magnetic field |l8j : nonlinear fiber lines with phase-sensitive amplification and 
mean-field models of degenerate optical parametric oscillators under continuous-wave pumping Il9l : pulsed optical 
parametric oscillators with spectral filtering and lasers with intracavity parametric amplification 20]. In most cases 
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the models considered in literature are either purely diffusive (c = oo) 0, 0, 0, or purely dispersive (c = 0) 
HI (HI E> EJ H> El- However there are situations where it is crucial that c be finite but nonzero. One physical 
phenomenon to which both diffusion and dispersion make essential contributions, is Faraday resonance in strongly 
dissipative media such as viscous fluids and granular materials plj . Another situation where < c < oo, corresponds 
to nondegenerate optical parametric oscillators |22| and fiber-optic telecommunication links; in these contexts, the 
term —icip xx represents spectral filtering in the spatial or temporal domain, respectively. In fact the applicability 
of Eq.@ with nonzero c and g may happen to be even wider; it is commonly held that this equation provides a 
phenomenological description to a broad range of pattern- forming systems |23( . 

To study the solitary wave phenomenology introduced by taking the (—icip xx )-teTm into account, we will use the 
diffusion- free limit (c = 0) as a starting point, and perform the numerical continuation of analytical solutions available 
in that case, to nonzero c. The stability of solutions obtained in this way will also be studied numerically. We will 
show that "switching on" the diffusion gives rise to a sequence of stable multihumped pulses occurring in the vicinity 
of a certain particular value of the diffusion coefficient, cn m = cu m (/i, 7). The closer c is to cu m , the greater is the 
number of multihump solutions coexisting at this c. The solutions with more than 5 or 6 humps describe a flat 
"plateau" (where tp is equal to the stable flat nonzero solution) sandwiched between two fronts of opposite polarity. 

The paper is organised as follows. Section [H] deals mainly with spatially homogeneous solutions. In particular, 
we show that there is a stable flat nonzero solution for sufficiently large c, and this uniform solution can serve as 
a background for solitary waves. In section UTTl we use perturbation- type arguments to demonstrate that both ip + 
and ip- solitons are continuable in c and in section IIVI we construct the solutions with nonzero c in the adiabatic 
approximation. The actual continuation is carried out numerically in section where we also examine the stability 
of the continued solutions. Some additional insight into the bifurcation of stationary pulses is gained in section I*VT1 
Finally, section TVIII summarises results of this work. 

We close this introduction by mentioning a recent publication [24j which was devoted to the study of the single- 
humped solution of Eq.@, by means of an averaging technique and direct numerical simulations. Neither the multi- 
stability of pulses nor the pulse- front transitions were dealt with in Ref . Q ■ 



II. EXISTENCE AND STABILITY INEQUALITIES 

The transformation ip(x,t) — > e~ luJt tfj(x,t) casts Eq.(j2J in a 'standard' autonomous form 

iipt + (1 - ic)ipxx + 2 \4>\ 2 ip - ip = hip* - ijip, (4) 

where we set g = and rescaled t so that u> = 1. This is the representation that we are going to work with in what 
follows. 

The equation 10} has three time-independent spatially uniform solutions, or 'flat backgrounds', for short. (We do 
not distinguish between solutions different only in the overall sign here.) One flat solution is i/'flat = 0; it will be 
central for this work where we are focussing on solutions decaying to zero at infinities. The other two flat solutions 
are given by 

^flat = *i 0) = (A ± /V2) e-* e ± , (5) 

where 

A± = \jl±^Jh 2 -7^, 26± = arccos(± v / l-7 2 //i 2 ). (6) 
Note that these solutions do not depend on c. 



A. Stability of spatially uniform solutions 



One can easily check that the zero solution is stable as long as h < yl 4- 7 2 , irrespectively of the value of c > 0. 
The analysis of the flat nonzero solutions is somewhat more laborious. 

Linearising Eq.Q) about Vflat — and assuming a perturbation Sip oc exp [i(uit — kx)] yields the dispersion 

relation 



{ck 1 + 7) ± WZ , (7) 



3 



where 

Z=(l-2A 2 + k 2 ) 2 + A 2 (A 2 -2)-h 2 . (8) 

(Here A stands for A + or A_, depending on which solution we are linearizing about.) The flat solution is stable iff 
Ke(ioj) < 0, i.e., when Z > —(ck 2 + j) 2 , for all real k. The latter condition amounts to an inequality 

(l + c 2 ) s 2 + 2 (l - 2A 2 + 7 c) s + AA 2 {A 2 - l) > 0, (9) 

where s stands for k 2 . 

Let, first, ip — . Since A 2 _ < 1, the inequality (O does not hold for k = 0; hence the 'low' background is 
always unstable. 

Let now tp = the 'high' background. The quadratic expression in @ will be positive for all s > if either its 
two roots si and S2 are both real negative or complex. (Note that we cannot have two real roots of opposite signs as 

the quadratic's constant term, AA 2 (A 2 — l), is always positive for ip — .) Whether the roots are real or complex 
is determined by the discriminant of the quadratic @ which can be written as 

V=[ 1 2 -AA\{A 2 + -mc-c + ){c-c_). 

Here we have introduced 

{2A\ - 1) 7 ± 2A + Ma 2 + -1){1 + ^) 
C± = 7^-4^(^-1) • (10) 

We need to consider two cases. Assume, first, that 4A 2 ^(A 2 h — 1) < j 2 . In this case the discriminant is negative 
(and hence, the roots s\ y 2 are complex) provided c lies in the interval c_ < c < c + . On the other hand, the roots s\2 
are real negative in this case if T> is > and the coefficient in front of the middle term in © is positive: 

c> c = (2A 2 + - 1) / 7 . (11) 

Since in the case at hand we have < c_ < c < c + , the union of the above two 'stable' regions, c_ < c < c + and 
c > c+, with the endpoints included, is simply c > c_. 

The second case is defined by inequality 4yl^_(A^ — 1) > -f 2 . Here we have c + < < c_ < c and the quadratic 
@ cannot have negative real roots as the region l|llfl does not overlap with the region where T> > 0. However, it can 
have complex roots — provided c > c_ . 

Thus we arrive at a simple stability criterion for the flat nonzero solution , valid for all h and 7: 

c>c_(/i,7), (12) 

with c_ as in (|10f> . Note that the stability threshold c-(h,j) is always strictly positive. This implies, in particular, 
that the solution is always unstable in the case c = 0. (This fact has already been mentioned in the literature 
|2^j.) Therefore the 'focusing', or 'attractive', damped-driven Schrodinger equation (Eq.Q with c = 0) does not have 
stable flat backgrounds except the trivial one, -0 = 0. The analysis of localised solutions over flat nonzero backgrounds 
becomes meaningful only within the full Ginzburg-Landau equation, i.e. Eq.Q with positive c. 



B. The flat solutions as backgrounds to solitary waves 

No less important is the question of when a flat solution can serve as an asymptotic value to a localised waveform 
— in other words, when is spatial decay to the flat solution possible. To find the corresponding criterion, we set uj = 
in equation (JJJ. This results in a quadratic equation 

(1 + c 2 ) ,s 2 - 2 (2A 2 - 1 - 7 c) s + AA 2 {A 2 - l) - 0, (13) 

where s — k 2 . The spatial decay to a flat background is possible unless both roots of Ea. i|13fl . si and S2, are positive. 

In the case of the ^ background, the discriminant of equation Ijl3(l is positive while the constant term is negative, 

whence s± > and S2 < 0. Consequently, the decay to may occur for all values of h, 7 and c. (This fact is of 
little importance, however, since we have just shown that this background is always unstable.) 
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In the case of the ^ + solution, we may have two positive roots — provided V > and c < cq, with cq as in 
Following the steps in the previous subsection, one can readily show that this situation occurs only if < c < c_, 
with c_ as in l|l(Jfl . Thus the stable constant solutions, defined by inequality (|12[) . can always serve as backgrounds 
to fronts and pulses. We will come across localised solutions over nonvanishing backgrounds in section Ivl below. 

Finally, to examine the case of the zero background, we set both A = and u> = in equation J7J). This yields 

(1 + c 2 )s 2 + 2( 7 c + l)s + 7 2 - h 2 = 0, (14) 

where, as before, s = k 2 . Since the middle term in l|14|) has a positive coefficient, we have si + 82 < 0, which 
means that the situation where both si and S2 are positive, can never occur in this case. Thus the decay to the zero 
background is possible for all h, 7 and c. 



C. The threshold driving strength for pulses 



Our last result in this section concerns the range of driving strengths which can support localised solutions in the 
presence of dissipation. Multiplying Eq.Q by ip* and subtracting the complex conjugate of the resulting equation 
gives what would be a local conservation law of the number of particles if c, h and 7 were equal to zero: 



t d t \^\ 2 + ^ x r-r^) x 
= h [(v/) 2 - - 2z 7 + ic (4> xx r + r xx i>) ■ 

Assuming that ip and ip x — ► as \x\ — > 00 and integrating 1151) over the real line, we get 

- ( \^\ 2 dx = 2 I |Vf[/isin(2 X )-7]^-2c f \^j x \ 2 dx, 



(15) 



(16) 



where we have denoted — x(x, t) the phase of the complex field ip: ip — \ip\e~ tx . Let h < 7 (and remember that c > 0). 
In this case it follows from (|16|l that the time derivative of / |i/'| 2 (ia; remains negative at all times. Hence as t — > 00, 
ifj{x,t) — > for all x. No stationary, time-periodic, quasiperiodic or chaotic solutions, decaying to zero as |x| — * 00, 
can arise if h < 7. We will make use of this simple criterion in what follows. 



III. CONTINUABILITY OF THE TWO NLS SOLITONS 



In the limit c = 0, Eq.Q becomes the parametricall y d riven damped nonlinear Schrodinger (NLS) equation, which 
has exact time-independent solitary wave solutions 0,^3 °f the form 

ip±{x,t) = A±e~ ie± sech(A±x) , (17) 

with A+ and 0± as in ffijl. The solution xp- is unstable for all h and 7, while "0+ is stable in a certain parameter 
region y|. 

The purpose of this section is to show that the solitary pulse solutions ip± of the NLS equation persist for c ^ 0. 
We restrict ourselves to stationary solutions (ipt = 0). Writing ip — (f>e~ lG with 9 a constant phase to be chosen later, 
Eq.Q) becomes an ordinary differential equation for <fi: 

(1 - ic) <\> xx + 2 |0| 2 4> - (1 - ii) <t> = h<j)*e 2m . (18) 

To find out whether solutions available at c = can be continued to nonzero c, we expand (f> in power series 
4> = 4>o + c 4>i + c 2 02 + • ■ ■) substitute into Eq. l(T%|) and match like powers of c. It is convenient to choose to be 6+ 
in the case of the soliton and 0_ for ip_; this choice makes (f> real. Matching terms linear in c and decomposing 
0i into its real and imaginary part (tfii = u + iv), yields an equation 

lJ:) = (_D, (19) 



where the primes stand for the derivatives in x and the operators L± are defined by 

/-a 2 + i-6^ 2 + / l cos(2e±) 27 . 

^±-1 -a 2 + i-2^ 2 -/ l cos(2e±) J ■ 
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(The subscripts + and — indicate whether we are checking the continuability of ip+ or ip—.) 

According to Fredholm's alternative, Eci. (|19|) has bounded solutions if and only if its right-hand side is orthogonal 
to the kernel of the adjoint operator L±. Transforming to £ = A + x and £ = A_x in the case of the ip + and ip_ 
solution, respectively, the equation for the zero modes spanning ker L\. becomes 



Li \ff 
27 L a - e± J \ g 



0, (21) 



where / = /(£), 3 = ff(0> 



e ± = 2 72 = ±2- 



4± 1 ± ^2 - 7 2 

and we have introduced the standard Poschl- Teller operators with familiar spectral properties: 

L Q = -a| + 1 - 2sech 2 £, L x = -9| + 1 - 6sech 2 £. 

The only discrete eigenvalue of Lq is Eq = 0, while the continuous spectrum occupies the semiaxis E > 1. Therefore, 
the operator (Lo — e±) is invertible as long as e± < 1 and e± 7^ 0. Assuming that the zero background is stable (i.e., 
assuming that h 2 < 1 +7 2 ), the quantity e+ is indeed less than one but greater than zero. On the other hand, the e_ 
is negative. (An exception is the point h = 7, where we have e± — 0.) Thus we conclude that (Lo — e±) is invertible 
except in the special case h = 7. (However, even in this special case the operator Lq is invertible on odd functions 
because the eigenfunction sech£ associated with the zero eigenvalue is even.) Consequently, /(£) cannot be equal to 
zero for h 7^ 7 — otherwise the bottom equation in l|21|l would imply that g(£) = 0, too. Fortunately, the operator L\ 
does have a zero eigenvalue, and therefore /(£) can be a nonzero multiple of the corresponding eigenfunction (which 
is tanh£ sech£.) 

Thus in the case h 7^ 7 the kernel of is spanned by just one zero mode, namely 



fy\ ( tanh £ sech £ 

91 J ~ \-2j(L - e±)~ 1 (tanh£sech£) 



(22) 



On the other hand, when h = 7, the dimension of the kernel space is two. Firstly, the zero mode persists for 
e± = as the operator Lq 1 is defined on odd functions. Secondly, there is another zero mode given by 



/ a W 

g 2 I \ sech£ 



(23) 



It is obvious that the vector-function (|22|l . whose both components are odd functions of £, is orthogonal to the 
right-hand side of Ea. (|19|l — which is even. Hence Ea. l|19|) is solvable for h ^ 7. It is also easy to check that the 
mode (12. '{t is not orthogonal to the right-hand side of (|19f) and so the solvability condition is not satisfied for h = 7. 

Thus, having started from the two nonlinear Schrodinger solitons, we constructed two one-pulse solutions of the 
Ginzburg-Landau equation to within the first order in the small parameter c: ip = ip± + ce _ * e± (u + iv) + .... 
Consequently, we expect to be able to continue the Schrodinger solitons ip± into the region c ^ (provided that 
h 7^ 7.) This expectation is born out by results displayed in section Ivl below, which present an outcome of the 
numerical continuation of ip± in c. 



IV. ADIABATIC APPROXIMATION 



Before attempting the full-scale numerical continuation, it is instructive to construct approximate solutions. Our 
approximation will be valid for small c and exploit the fact that when c = 0, the stationary pulse- like solutions of 
equation J2J have the form 

ip = asech{ax)e~ ze , (24) 

where a = A± and 9 = Q± are constants defined by Eqs.©. For c small nonzero, approximate solutions can be 
sought for by assuming that ip retains the form l|24|l . but a and 9 become slowly varying functions of t. 
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To obtain an expression for a (the overdot stands for the derivative in t), we substitute the ansatz (|24(l into 1(15(1 . 
Integrating over x and using the boundary conditions ip x — > as \x\ — > oo produces an evolution equation for the 
pulse's amplitude: 

a = 2a(h sin 26- 7- ca 2 ) , (25) 

where c = c/3. An equation for the pulse's phase arises by multiplying by ip* , adding with its complex conjugate, 
substituting (|2fjl for ip and integrating over x: 

9 = hcos29+l- a 2 . (26) 

Fixed points of the system (|25|l and ((261) correspond to stationary solutions of (@J. These can be easily found 
explicitly: Eliminating 9 from the stationary system 



produces a quadratic equation 



with roots 



h sin(20) = ca 2 + 7, /i cos(20) = a 2 - 1, (27) 
(1 + c" 2 ) a 4 + 2 (7c - 1) a 2 + 1 + 7 2 - h 2 = (28) 
2 _ (1 - 7 5) ± V(h 2 - 1)5 2 - 2 7 5+ fe^ 2 " .... 

a± r+i 2 • (29) 

The corresponding 9± are defined by their sine and cosine in 127(1 . 

It is not difficult to determine when the roots 1)29(1 are real and positive. We assume that the zero background 
solution of Eq.QJ is stable, i.e. the constant term in 1(28(1 is positive. Hence if the roots are real, they are of the same 
sign, and this sign is opposite to that of the middle term in ((28(1 . Therefore we have two positive roots provided the 
discriminant 

V = (h 2 - 1) c 2 - 2 7 5 + h 2 - 7 2 = (h 2 - l) (5 - ci)(c - c 2 ) (30) 
is nonnegative, and, at the same time, the inequality 7c < 1 holds true. In l(3U(l we have introduced 



7 ± Vl+ 7 2-fo2 

c i- 2 - — w^i — ' (31) 

where the + corresponds to c\ and — to c.%. 

It is straightforward to verify that for small h, h 2 < 1, we have c\ < < 62 < while for larger h, h 2 > 1, we 

have < 6% < — < cy, (Here we have made use of the inequality h > 7; as proved in section fll CI no stationary 
localised solutions exist for h < 7.) Therefore, the region of c values where the roots of 128(1 are positive, is given by 
the inequality c < c^Qi, 7) — for all h. 

Thus we conclude that the adiabatic equations have two stationary points for c < 3c~2(/i, 7), and none for c > 3c~2. 
We complete the adiabatic analysis by classifying their stability and bifurcation. 

Linearising equations ((26(1 about the stationary points and assuming small perturbations of the form 5a = Cie 2Xt 
and 89 = C2e 2At , yields a characteristic equation 

A 2 + (7 + ca 2 )A + 2ha 2 cos(2(9) - h sin(20) [h sin(20) - 7 - ca 2 ] = 0. 

Since the coefficient in front of the middle term in this equation is positive, either its two roots are complex with 
negative real parts, or we have two real roots, of which one is negative. The case when the second root is positive 
(and hence the fixed point is unstable) occurs if the constant term is negative. Conversely, if the constant term is 
nonnegative: 

2ha 2 cos(26>) - h sin(20) [h sin(20) - 7 - ca 2 } > 0, (32) 

the fixed point is stable. Simplifying Eg. 1(32(1 . we get the stability condition in the form 

2 ^ 1 - 7c 
1 + c? 
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FIG. 1: The bifurcation diagram displaying the Sobolev norm versus the diffusion coefficient, for the solitary pulse obtained by the 
continuation of the stable soliton in c. The continuation starts from a point on the (c = 0)-axis, marked by an open circle. The 
arrows indicate the directions of continuation and are only added for reference purposes. Solutions at points marked by the black 
dots (a), (b) and (c) are shown in Figure [5] The solid lines correspond to stable and dashed curves to unstable branches. 



Comparing this to the expressions for a± in (|29H . we conclude that for all c, h and 7, the points (a + , 8 + ) and (a_, 0_) 
are a stable node and saddle, respectively. 

The two fixed points come into being through a saddle-node bifurcation which occurs as the diffusion coefficient c 
is decreased past c = 3c2 for the fixed h and 7 or, alternatively, as the driving amplitude h is increased for the fixed 
dissipation coefficients c and 7. One can easily find the bifurcation value of h, at which two complex roots of the 
quadratic equation (|28|l converge on the positive real axis. (Here we are assuming that c is smaller than -.) Writing 
the discriminant l|30|) as 



V = {c 2 + 1) 



(5 + 7) 



- 1 



the threshold driving strength is given by 



had 



c + 7 



(33) 



where the subscript "ad" serves to remind that Eq. l|33fl was obtained in the adiabatic approximation. It is important 
to emphasise that this formula is valid only for small c. Equivalently, the expression for the turning point c 2 is 
valid only for h close to 7. 



V. NUMERICAL CONTINUATION AND STABILITY ANALYSIS 



In this section we describe the bifurcation diagram obtained by the numerical continuation of the solitons ip± in 
the parameter c. The diagram is presented in Fig^ and displays the Sobolev norm of the solution, 



as a function of c. 



A. The method 



For stationary solutions, ijj = ^(x), equation Q reduces to an ordinary differential equation 

(1 - ic)ip xx + 2 |-0| 2 $ - (1 - ii)1> = h^*- (34) 

For the numerical continuation of solutions to Ea . (|34(l we utilised the AUT097 software package j2||. The infinite 
line was approximated by a finite interval (—L,L), with L — 100, and the boundary conditions ip(±L) = 0. The 
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tolerance of the computation (the grid norm of the difference between the right and left-hand side of 134| l) was set to 
be 1(T 7 . 

The solid and broken curves in FigQ represent stable and unstable branches, respectively. The stability was 
examined by linearising Eq.Q about the corresponding stationary solution. Choosing the small perturbation in the 
form dtp(x, t) = [u(x) + iv(x)] e At , with real u and v, we arrive at an eigenvalue problem 

^(:)= Aj (:)< (35) 

where 



/ -dl + 1 + h - 6TZ 2 - 2X 2 -cdl + 7 - 4^2" \ 

n ~ \ cd 2 x - 7 - AKX -81 + 1 - h - 2K 2 - 6X 2 J ' 

and 




The 1Z and X are the real and imaginary parts of the solution: ip = TZ + iX. We solved the eigenvalue problem by 
expanding u and v over 500 Fourier modes in the interval (—100, 100). 



B. Continuation and stability 



The continuation in c was performed for fixed values of h and 7. We selected 7 = 0.5 and h = 0.8; for these h and 
7 the nonlinear Schrodinger soliton tp + is stable |8j. Before proceeding to the ip+ soliton, however, we briefly deal 
with the ^--case. In agreement with predictions of section ITTT1 the soliton ip- was found to persist both for c < 
and c > 0. We were in fact able to continue it indefinitely without encountering any obstacles in both directions. As 
c — > 00, the width of the localised solution of i|3~4*)l grows and the solution tends to ^^(x/y/c), where </> 7i h,(X) stands 
for a pulse-shaped solution of equation 

-icbxx + 2\<p\ 2 (P- {l-i 1 )<t> = h<j>*. (36) 

In a similar way, as c — > — 00, the solution tends to <f>* h {x/\f^c). When continued to c > 0, the pulse remained 
unstable for all c, with a single positive eigenvalue in its spectrum. (As for the negative-c region, all solutions there 
are a priori unstable against continuous spectrum perturbations with arbitrarily large ReA.) 

The continuation of the ip+ soliton proved to be more rewarding from the stability viewpoint. The corresponding 
bifurcation diagram is displayed in Fig^ As with the soliton the ip+ persists both for c > and c < — in 
agreement with section IIIII Continuing into the c > region, we have found that the solution gradually changes 
its shape, with the hump in the imaginary part splitting into two (see FigEJa)). At c = 0.7 the branch turns back 
(FigHJ. Note that the tur ning point occurs for c much smaller than 3c2 = 1.04, the saddle- node bifurcation point 
predicted by the adiabatic analysis. What is more important, the ip+ solution turns not into the ip- (as the adiabatic 
approach predicted) but into some other solution which has two well-separated humps in the imaginary part. The 
discrepancy is not surprising as the adiabatic approximation is valid only for small c, where the shape of the solution 
is still well reproduced by the single- hump constant-phase trial function (|24|l . In the entire stretch between c = and 
the turning point, the solution remains stable. 

Continuing this branch further, additional humps are added as the solution passes a sequence of turning points. 
Each pass of a turning point results in the creation of a new hump in the middle of the solitary wave. As we move 
along the branch, the distance between successive turning points (the difference between the corresponding values of 
c) becomes smaller and the new humps come with smaller amplitudes. As a result, a long plateau is formed which 
keeps on expanding as we continue the branch (FigP|c)). The broadening plateau accounts for the vertical segment 
of the curve in Fig^ If the bifurcation parameter c is seen as a function of the Sobolev norm \ \ip\\ (that is, if we turn 
Figfflby 90°), the curve c(||^||) has the form of a decaying oscillation. 

The stability of the solution alternates at each successive turning point. These changes are due to a single real 
eigenvalue which moves back and forth through the origin on the real line. (At the turning points the eigenvalue is 
right at the origin, of course.) The lengths of the incursions this eigenvalue makes into the positive and negative real 
lines decrease with each new turning point until the eigenvalue becomes indistinguishable from zero. Therefore the 
branch becomes (neutrally) stable sufficiently 'high up' in \ \ip\ \ in Fig^(i.e. for sufficiently long plateaus.) 
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FIG. 2: (a),(b),(c): Solutions at the corresponding points in Fig0 Solid line: real part; dashed line: imaginary part. 

Despite the fact that negative values of the diffusivity c are not physically meaningful, we did continue to c < — 
in the hope that the resulting branch would reach a turning point and then return to the positive c-semiaxis (which 
it did indeed do). The resulting branch is also shown in Fig^ As we move into the region c < 0, the ip+ solution 

gradually develops into a three-humped state and when the curve returns to c = 0, we have a complex of the | \ 

type, with a large separation between the individual solitons. (This complex of three solitons of the parametrically 
driven nonlinear Schrodinger equation was previously found in Ref.jjJ.) Crossing into the c > region, the central ip+ 
soliton in the complex transforms as if it did not have the ?/>_ solitons attached to its flanks. As a result, the c > 0- 
portion of the corresponding ||^(c)|| curve has virtually the same shape as the curve resulting from the continuation 

of ip+ to the region c > 0; the only difference is that the curve emanating out of ip^ | ) is shifted upwards relative 

to the curve emanating out of tp+. Similarly to the continuation of the continuation of ^( | ) goes via a series 

of turning points, with each pass of the turning point resulting in the creation of another hump in the middle of the 
central region which becomes a long plateau. The lateral ip_ solitons are not affected by this process. The linearised 
spectrum is the union of the spectrum of the long pulse described in the previous paragraph and spectra of two -0_ 
solitons. In particular, it includes two positive real eigenvalues contributed by the V'-'s and so the entire branch 
resulting from the continuation of ipr | ) is unstable. We disregard it in what follows. 

The plateau arising in the final stage of continuation of the two branches shown in Fig^ is nothing but an interval 
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on the x-axis where ip equals $ L , the flat nonzero solution given by Eq.lJSJ. The corresponding value of c, cy im = 0.54, 

falls within the region c > c_ where the background is stable. Here c_(7,/i) is given by Eci. (|10|) ; in particular, 
c_ (0.5, 0.8) = 0.224. In the language of phase transitions, the long pulse shown in Fig|2fc) can be seen as a "bubble" 
of one stable phase in another one. 



The long pulse can also be interpreted as a bound state of two fronts interpolating between different stable back- 
grounds, ip — and ip — . This intuitively-appealing idea can be put on the quantitative footing by considering 

the spatial decay of perturbations to the flat background 'J^' . 

Indeed, consider the turning point separating the second stable branch from the first unstable branch in Fig^ 
(This is a point with coordinates c = 0.37 and ||-0|| = 2.48. Note that we are only considering branches obtained by 
the continuation of the soliton tp+ to positive c. Branches obtained by the continuation of ip- to negative c first, are 
disregarded here.) It is at this point that the modulus squared of ip{x) becomes double-humped; before that, that 
is on the branch that leads to this point (the first unstable branch), the function |?/;(a;)| 2 remained single- humped 
despite the double-humped imaginary part. Letting h = 0.8, 7 = 0.5 and c = 0.37 we check that 4A 2 h (A 2 h — 1) > 
and hence, according to section lll Al the quadratic 10, 1|13|) has two complex roots, k 2 = s r +isi and (fc 2 )*. Therefore 
the wavenumber k is complex as well: k — k r + i/cj. Solving equation i|13|) for s we obtain, subsequently, 



(We have chosen positive values for k r and ki here.) 

For c near the turning point in question, the plateau has not yet formed between the two humps and so they can 
be crudely thought of as two overlapping NLS-like solitons with oscillations on their adjacent (i.e. partner-facing) 
tails. The bound state arises when one solitary wave is trapped in the potential well formed by the oscillatory tail of 
its partner. Making use of the potential of interaction of two attractive NLS solitons j^], U e g oc exp(— kiz) cos(k r z) 
(where z stands for the distance between the two humps), and taking into account that kijk r <C 1, we obtain a rough 
estimate for the separation: z — n/k r . This formula could also be derived from purely qualitative considerations; all 
one needs to notice is that the absolute value of the field ip obtained from the linear superposition of two in-phase 
solitons is maximised (and so the energy is minimised) by placing one soliton at the first maximum of \ip\ on the tail 
of its partner. Substituting k r from (|37() . the approximate formula gives z = 2.31. This is in qualitative agreement 
with the numerically found value z = 2.77. 

Moving further along the (second stable) branch, the plateau appears and the two humps can no longer be approxi- 
mated by the NLS-like solitons. This makes the above estimate invalid. The solution can still be regarded as a bound 
state of two fronts but this time, in order to calculate the characteristic width of the pulse one would need to know 
the full profile of the front. 

We conjecture that for a given h and 7, a free-standing stationary front exists just for a single value of c, namely 
c = ci; m . (We are planning to verify this conjecture in future publications.) On the other hand, stable and unstable 
bound states of fronts exist in a finite interval of c values containing c = cn m as an internal point. It is fitting to note 
here that similar pulse-to-front transformations occur also in the other system featuring the subcritical bifurcation, 
viz. the cubic-quintic Ginzburg-Landau equation with internal gain ps| . 



As we mentioned in the previous section, there is a certain discrepancy between the adiabatic analysis and numerical 
continuation. Numerically, the ip- soliton was found to be continuable all the way to c = +00 whereas the adiabatic 
approach predicted the existence of a turning point at c = 3c2, where the ip- should have merged with the tp+ branch. 
As for the ip + solution, we found that it turns into a pulse with the double-humped imaginary part (and not into the 
ip- branch as suggested by the adiabatic approximation.) In order to shed some light on the possible source of the 
discrepancy we performed the numerical continuation of the pulse ip- in h, for several fixed values of c. Here, by the 
ip- we mean the pulse solution which results from the continuation of the Schrodinger ip- soliton to positive c, for 
some fixed large value of h (in our case for h = 0.8.) Having obtained this starting-point solution for several values 
of c, we then continued it in h, from h = 0.8 to smaller h. The stability of the arising solutions was examined by 
computing eigenvalues of the operator (|35(l at sample values of h. 



C. The bound state interpretation 




(37) 



VI. THE h VS c DIAGRAM 



11 




FIG. 3: (a) The Sobolev norm of the solitary wave solution as a function of the driving strength h. The solid and dashed lines 
indicate stable and unstable branches, respectively, (b) The pulse existence region on the (c, /t)-plane. Two pulse solutions, ip+ and 
tj)-, are born as h exceeds the value h CI (c) depicted by the lower solid line. The upper solid curve gives the upper boundary of the 
T/>+-pulse's existence domain, /12(c). Also shown is the adiabatic approximation to the saddle-node bifurcation curve, Ea.l)33[l (dashed 
line). 



In each case considered, the ip- branch was found to turn into the ip+ solution as h reached the threshold value 
h cr = h cr {-^.c). (That is, for h > h cr there are two branches of solutions whereas for h < h CT , there is none; see 
FigEfa)). The entire ip- branch is unstable; the single positive real eigenvalue moves to the negative semiaxis as the 
branch is continued past the turning point. Continuing the arising ip + branch to larger h, we reach another turning 
point at h = h,2, where the ip+ solution transforms into a pulse with the double- humped imaginary part. The values 
h CI and h,2 are shown in Fig[3Jb) , as functions of c (for the hxed 7 = 0.5). As c — > 00, the difference h% — h CT decreases 
but remains nonzero. We verified this by computing h cr and h 2 for equation i|36|l which pertains to c = 00. In the 
same plot we display the function which gives the adiabatic approximation to the curve h CI (c). (Note that for 
small c, there is a good agreement between numerical and approximate values but as c grows, the two curves diverge.) 

Continuing the ip+ branch past the second turning point, the solution adds another hump in the middle of the 
pulse, turns back again, adds another one, and so on. (See Fig|3fa)). A long plateau developes in the middle of the 
pulse, just like when it was continued in c. Similarly to the c-continuation, the turning points on the ft,-axis separate 
regions of stability from regions of instability, with the instability being caused by a single positive real eigenvalue 
(which moves through A = at the turning points.) 

From figure EJ D ) it is clear why the saddle-node bifurcation point where the ip + and ip- solutions would merge, 
did not appear in Fig^ The reason is that FigQ] was plotted for a relatively large value of h (h = 0.8) whereas 
according to Fig|3{b), a horisontal line h = const with h > 0.660 can have no intersections with the saddle-node curve 
h CI (c). (Here all the numbers are for 7 = 0.5.) The same Fig|3Jb) explains what semed to be a discrepancy between 
the adiabatic analysis and the numerical continuation of the soliton ip- in c. The numerical result that seemed to 
contradict the adiabatics was that for h = 0.8, i/j- could be continued without bounds. It is now obvious from Fig|3Jb) 
that the unbounded continuation is only possible for h greater than 0.660. Continuing the ip- soliton to positive c for 
h smaller than 0.660, the branch turns back (already as the pulse) after hitting the lower solid curve in Fig|3^b). 
Therefore, the pattern arising for h close to 7 actually is in qualitative agreement with the adiabatic analysis, which 
was expected to be valid precisely for small c or, equivalently, for small h — 7 differences. 
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VII. CONCLUDING REMARKS 

In this paper we studied a cubic complex Ginzburg-Landau equation in which linear losses and diffusion are com- 
pensated by the linear parametric drive. The nonlinear term in the equation was taken to be purely conservative. 

There are three stationary homogeneous solutions to equation (0J, and we have shown that the ip = solution is 
stable if h < yl + 7 2 , as long as c > 0. This stability condition coincides with the corresponding condition for the 
nonlinear Schrodinger case (i.e. for c = 0). On the other hand, the stability properties of the nonzero homogeneous 
solutions are not the same as in the c = case. Indeed, unlike for c = 0, there is a stable flat nonzero solution 
if) = \fj^ for sufficiently large diffusion coefficients, c > c^(h,"f). 

Having established the persistence of the NLS solitons -0+ and ip- for small nonzero c, we continued them in 
c numerically. The continuation of the soliton ip + yields a sequence of pulse-like solutions, separated by turning 
points, with increasing number of humps. The stability of these solutions changes at each turning point, so that 
stable multihump solutions coexist with unstable ones. It is fitting to note here that the multistability of multipulse 
solutions is not observed in the Schrodinger limit where only the two-soliton complex was found to be stable |9|. As 
c — > Qim. where cn m — ci; m (/i, 7), the solution takes the form of a long plateau (an interval of the stable background 

sandwiched between two fronts. 
We also performed the continuation in h, for the fixed c. For each c > 0, two localised solutions are born in a 
saddle-node bifurcation as h exceeds a threshold value. (We obtained an analytic formula for the threshold in the 
adiabatic approximation; it is in agreement with the numerical results for small c.) The subsequent continuation 
gives rise to a sequence of coexisting stable multihump solutions culminating in a bound state of two widely separated 
fronts. 

Solitary pulses in the form of long shelves (plateaus) can allow easy experimental observation in physical systems 
described by our model. In particular, they may be employed as a natural basis for the non-return-to-zero (NRZ) 
format of the data transmission in optical telecommunications. In the NRZ format, the 1 and bits are coded, 
respectively, by sending or withholding the signal within a standard time slot. 

A string of several l's looks as a long uniform pulse of an essentially arbitrary length. The stability of such pulses 
is crucial to maintain the fixed shape of the long array of l's, and to prevent the inter-symbol interference, i.e., the 
blurring of empty intervals between such strings, which represent (strings of) O'sfsee, e.g. Ref. j2|| and references 
therein.) In the case of lasers, which can also be described by the present model [2(J], the possibility of the generation 
of stable long pulses of arbitrary duration, i.e., an effective tunability of the output, is an essential advantage too. 
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